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AUTOMATIC STABILIZING SYSTEM 

EXPERIMENTAL AND CALCULATED DYNAMIC STABILITY 
CHARACTERISTICS FOR SPEEDS UP TO 200 MPH 
By William C. Schneider 

SUMMARY 


Tests of a toved airplane model have been made in the langley 
19 “foot pressure tunnel in conjunction with the development of a flutter 
testing technique. The model was equipped with a suitable autopilot to 
keep the model flying strai^t and level in the tunnel while restrained 
only in drag. Flights were made through the speed range of the tunnel 
up to a maximum speed of 200 mph at atmospheric pressure with both rigid 
and flexible wings. 

The dynamic stability of the model-autopllot-towline combination 
was investigated analytically and sui attenq)t was made to correlate the 
results with experiment. Only qualitative comparison with flight 
iii5)ressions was possible due to a randcm tunnel airstream distiirbance 
of unknown origin. 


INTRODUCTION 


Tests in the langley 19-foot pressiore tiannel of a towed airplane 
model have recently been completed in conjunction with the development 
of a testing technique which allows body freedoms during flutter tests 
(ref. l) . The model was flown on the end of a tow cable attached to 
the nose of the model. Sxiitable autopilots were installed in the model 
for the purpose of keeping the model flying straxght and level in the 
tunnel. 
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In order to provide a statle test vehicle, the autopilot system 
had to perform adequately as a regulatory system which would keep the 
model flying within the confines of the tunnel, hut still would allow 
the model to behave largely as a free body during flutter tests . 

The model tested was a dynamically scaled fighter-type airplane 
with a wing of aspect ratio 5-45, 40° sweepback of the quarter-chord 
line, and a taper ratio of O.578. In order to simplify the initial 
experimental tests, .aeroelastic effects were avoided by using a set of 
rigid wing panels. An Interchangeable set of flexible panels was pro- 
vided for the wing flutter tests. 

Flights were made in the 19-foot pressure tunnel at atmospheric 
pressure up to a maximum airstream velocity of 200 nph. 

It is anticipated that the towed-model testing technique will be 
utilized to a greater extent than heretofore and that there will then 
be a need for information concerning the technique. This report is 
intended to demonstrate the feasibility of wind-tunnel testing by the 
towed-model technique at moderate subsonic speeds. Comparison of 
experimental lateral and longitudinal time histories of the model motion, 
for the model equipped with both the rigid wing panels and the flexible 
wing panels, is made with motions calculated on an analog computer. The 
lateral eqmtions of motion for a towed airplane have been developed in 
reference 2, suid the longitudinal equations of motion are developed in 
an appendix to the present paper. This comparison was made despite the 
knowledge that experimentally the model was excited by random disturb- 
ances of an unknown magnitude, so that the coo^parison was not orthodox. 

In euidition, lateral periods and dajjiping calculated for various auto- 
pilot conditions are included. 


SYMB0I5 


All forces and moments are referred to the stability system of axes 
which is described in figure 1. 

A,B,C,E,F,G coefficients of the lateral-stability equation 

H,I,J,K coefficients of the longitudinal-stability equation 



drag 

coefficient. 

D 

q'S 


lift 

coefficient. 

L 

q'S 
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■'in 


pitching-moment coefficient. 


M 


q'Sc 


Cy 


Cl/2 

D 

Db 


T>r 


% 




ix 


Iz 


yawing-moment coefficient. 


rolling-moment coefficient. 


N 


q'Sb 


q'Sb 


weight coefficient. 


W 

q‘S 


towline longitudinal-force coefficient. 


side-force coefficient, 

q'S 

Tt, 

towline vertical-force coefficient, 

q*S 


cycles to deunp to half anqplitude. 


drag force , lb 


Ti/2 


lateral differential operator, 

ds-b 


longitudinal differential operator, 

dsc 

moment of inertia about principal X-axis, sliog-ft^ 
moment of inertia about principal Y-axis, sliig-ft^ 
moment of Inertia about principal Z-axis, slug-ft^ 
mcanent of inertia about X-axis, slug-ft^ 
moment of inertia about Y-axis, slug-ft^ 
moment of inertia about Z-£ixis, sliig-ft^ 
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product of inertia^ - ^IZq " ^ ^ 


ix ^ 

radius of gyration atout X-axis^ — 

mb^ 

k/ 

ly i 

radius of gyration about Y-axis^ 


I 7 I 

radius of gyration about Z-axis, — ^ 

mb^ 

KxZ 

product-of -inertia factor, 

mb*^ 

CD 

elevator -position-control gearing ratio, radians/radian 


pitch-daraper gearing ratio, radians/radian/sec , 


roll-autopilot gearing ratio, radians/radian 

K^i,' 

• 

rudder-position-control gearing ratio, radians/radian 


yaw-damper gearing ratio, radians/radian/sec 

L 

lift, lb j 

L* 

rolling moment, ft-lb 

M 

pitching moment, ft-lb 1 

N 

yawing moment, ft-lb 

P 

period of oscillation, sec , 

S 

wing area, sq ft 

T 

towline tension, Ih 

Ti/2 

time required for oscillation to damp to one-half 
ajnplitude, sec 

V 

tunnel airstream velocity, ft/sec 

w 

model weight, lb 
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Y 

b 

c 

c 

I 

m 

P 

q 

q' 


r 


Sb 


side force, lb 

model wing span, ft 

model wing chord, ft 

model wing mean aerodynsunic chord, ft 

towline length, ft 

mass of model, sliigs 


rolling velocity. 



pitching velocity. 


dt 


tunnel dynamic pressure. 


i pV^, Ib/sq ft 


yawing velocity. 


dT|f 

dt 


nondimens ional time parameter based on span. 


Vt 

b 


Sq nondimens ional time parameter based on mean aerodynamic 

chord, — 
c 

t time, sec 

u nondimens ional longitudinal perturbation velocity 

x' distance along model center line from center of gravity 

to towline attachment point, positive forward of 
center of gravity, ft 

y lateral displacement of model center of gravity, positive 

to the right, ft 

z vertical displacement of model center of gravity from 

equilibriimi, positive down, ft 

z ' vertical distance from X-axis to towline attachment point 
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perturbation angle of attack, radians 
angle of attack, radians 

angle of attack in undistiarbed condition, radians 

rate of change of angle of attack, — 

dt 

angle of sideslip, radians 
fli^t-path angle, radians 

angle of aileron deflection perpendicular to aileron 
hinge line, positive to produce positive roll, radians 

angle of elevator deflection perpendicular to elevator 
hinge line, positive trailing edge down, radians 

angle of rudder deflection perpendicular to rudder 
hinge line, positive trailing edge left, radians 

inclination of principal longitudinal axis, radians 

inclination of principal longitudinal axis from fuselage 
center line, radians 

perturbation angle of pitch, radians 

angle of pitch in undistiorbed condition, radians 

angle between towline and model center line in xz-plane, 
radians 

rate of change of angle of pitch, — 

dt 

angle between towline and horizontal, radians 

relative density factor based on span, -S- 

pSb 

relative density factor based on chord, 

^ pSC 

angle between towline and X-axis, radians 
tunnel airstream density, slugs/cu ft 
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angle of roll, radians 
angle of yaw, radians 
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0 

Hr* 

n 


rate of change of angle of yaw. 


di 

dt 


angle between towline and model center line in xy -plane, 
radians 

total angle of pitch, Qq "*■ radians 


Aerodynamic coefficients and force coefficients with a subscript indicate 
the partial derivative of the quantity with respect to the indicated 
quantity nondimens ionali zed in the standard manner. For exanqjle. 


C 


^Ci 


\2V dt/ 


MDDEL 


The tests were made on a dynamically scaled model of a swept-wing 
fighter-type aircraft of current design. Pertinent model areas and 
dimensions are shown in figure 2 and are listed in table I. Table I 
also lists the scaled mass characteristics of the model. The f\iselage 
and tail surfaces were constructed of a fiber glass and plastic laminate, 
wrapped aroiaid a shell of foam plastic for rigidity. The resultant 
structure was li^t, strong, and essentially rigid, and provided the 
housing for the autopilot system and the instrximentation. 

Two interchangeable pairs of wing panels were provided for the 
model. One pair of panels was rigid and was constructed of alumi num 
sheets sandwiched between mahogany layers. The other pair of wing peinels 
(described in detail in ref. l) was constructed to have scaled stiffness, 
weight, and inertia characteristics and was, as a result, flexible. The 
bending and torsion strength of these flexible wing panels was concen- 
trated in an alxominum spar. The airfoil contour was formed by independent 
balsa-wood pods which were attached to the spar. Full-chord wing fences 
were used on the model to improve the flow over the wing at low speeds. 

The fences were located at 68 percent of the wing span. 

External stores were mounted below the wings as shown in figure 2. 
These stores provided a housing for gyroscopic flutter stoppers which 
provided mechanical damping to the wing at the operator's command. The 
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gyroscopes were air-driven and were mounted in gimbals normally free to 
move. A solenoid-operated air valve permitted the operator to restrain 
the gyroscope gimbals which provided a damping force which opposed the 
wing flutter motion. 


AUTOPILOT AND CONTROIS 


Several mechanisms were incorporated in the model to keep it flying 
straight and level in the tunnel. Each of the mechanisms was independent 
of the other components, but all of the mechanisms as a group will be 
referred to as the autopilot in this paper. Each of the mechanisms 
individually will be discussed as a component of the autopilot. 

Roll control .- To provide roll control, an electrically driven dis- 
placement gyroscope was mechanically linked to the ailerons. The linkage 
was such that any bank angle caused an aileron deflection which provided 
a restoring rolling moment. Mathematically, the roll con5)onent of the 
autopilot could be represented by the eqxiation 6^ = where is 

the gearing ratio is the angle of roll about the model longitudinal 
axis) . The value of the gearing ratio could be varied, but was fixed 
during any particular fli^t. A photograph of the gyroscope is shown in 
figure 5^ and the gyroscope and linkage is shown schematically in figure 4. 
A description of the gyroscope can be found in reference 5 and the dynamic 
characteristics are presented in table II. 

The gyroscope tended to drift away from the zero position owing to 
both friction in the bearings and aileron hinge moments. This gyroscopic 
drift produced an aileron deflection with a resultant model angle of bank. 
In order to compensate for the drift during flight, an electromagnetic 
precession coil was mounted so that a corrective torque could be applied 
to the inner gimbal of the gyroscope at the operator's command in such a 
manner as to cause the gyroscope to return to the zero position. 

The ailerons were plain flaps with a chord of 25 percent of the wing 
chord perpendicular to the quarter chord. The aileron span extended from 
0.155b/2 to 0.517b/2. 

Automatic position control . - The towline was linked to the elevator 
and rudder in such a manner that a control-surface deflection was produced 
whenever the model center line was not alined with the towline. The 
linkage was such that the deflection produced a corrective moment to 
realine the model center line with the towline (see fig. 4). Mathemati- 
cally, the deflection could be represented by the equations 6g = K0i9' 

and 6j. = K^ii' where Kq' and K^' are the gearing ratios and 9' 
and i|f' are the angles between the towline and the model center line in 
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the yz- and xy-planes. The values of Kgi and K^, could he varied, 
but were fixed during any flight. 

The elevator segment connected to this linkage had a chord of 
50 percent of the stabilizer chord and had a span extending from 7 to 
50 percent of the stabilizer span. Similarly, the rudder segment 
connected to this linkage had a constant chord of 2 . 7^1 inches parallel 
to the fuselage center line and a span extending from 0.27 to 50 percent 
of the vertical tail span. 

Pitch and yaw dampers . - The model aerodynamic danqping in pitch and 
damping in yaw were augmented by rate-gyro~servo imits . These units 
provided elevator or rudder deflections proportional to the pitching 
and yawing velocities, respectively, which caused moments which opposed 
the motion, thus providing artificial damping. Uie dampers can be repre- 
sented mathematically by the equations 5 g = Kg 9 and 5 ^. = where 

and are the gain constants, and 0 and ^ are the pitching 

and yawing velocities. The values of Kq and K^. could be varied 

during fli^t. Schematic representation of the rate-gyro — servo units 
and the linkages are shown in figure 4 - and a description can be found 
in reference 4 . The dynamic characteristics of these units are listed 
in table II. 

The pitch damper was linked to the outboard segment of the elevator 
(span extending from 50 percent to 100 percent of the horizontal stabi- 
lizer span) and had a chord of 50 percent of the stabilizer chord. 

The yaw damp er was linked to the upper segment of the rudder (span 
extending from 50 percent to 88 percent of the vertical stabilizer span) 
and had a chord of 2.741 inches, parallel to the fuselage center line. 

Manual vertical position control . - During airspeed changes , the 
model was kept at or near the tunnel center line (vertically) by varying 
the deflection of the inner segment of the elevator (the segment linked 
to the automatic position control) . A small electric motor was incor- 
porated in the automatic position control linkage in such a manner as to 
vary the length of a push rod, at the operator's command, thereby changing 
the elevator deflection. This manual control was also used to take off 
and land the model. 


TESTING EQUIPMENT 

The towUne was 58 feet long (6.8b) from the point of attachment at 
the nose to the restraining yoke. The yoke was in the entrance cone of 
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the tunnel and was restrained longitudinally by a drag cable attached to 
the timnel turning vanes. (See fig. 5 *) Vertical and lateral motions 
of the yoke were restrained by three cables which fixed the yoke on the 
tunnel center line. 

The towline consisted of 55 feet of l/l6-inch flexible aircraft 
cable plus a 5-foo'fc aluminum rod attached to a universal joint at the 
nose of the model. The al umi num rod was required to provide a moment 
at the universal joint to oppose the control-surface hinge moment. A 
bundle of electrical leads and a flexible plastic air supply tube was 
attached to the aircraft cable from the yoke to the model. The leads 
provided electrical power to the model and transmitted data from the 
instruments. The air supply tube provided air which was required for 
the flutter-stopper gyroscopes and for the rate-gyro — servo vinits. The 
total weight of the towline, lead.s, and air tube was about 5-5 pounds. 

A flat rubber padded platform provided a landing mat. (See fig. 5 -) 
A bubble canopy was installed in the tunnel floor to the rear of the test 
section to provide an observation point for the aileron trim operator. 

As a safety feature, a model hold-down wire was provided to keep the 
model on the landing mat before and after the flights. This wire con- 
sisted of l/l6-inch aircraft cable attached at the bottom of the fuselage 
just below the center of gravity. The cable was conducted out of the 
tunnel to an operator who provided the pull required to keep the model 
on the floor. During fli^t, the wire was kept slack. 

The model motions were recorded by three motion-picture cameras 
placed at the top, side, and rear of the tunnel test section. Instru- 
mentation in the model provided information as to the model accelerations 
and the flexible wing panel motions (see ref. l) and the control-surface 
deflections. These data were recorded on an oscillograph but are not 
presented. 

A timing device (which indicated on the camera film and on the 
oscillograph records) enabled the two sets of data to be correlated with 
respect to time. 

Data were recorded at 100 mph and at the maximum speed of each 
flight. 


TESTING PROCEDURE 


The model was placed on the landing mat in a position which put 
tension into the tow cable approximately equal to the model drag at 
take-off. All controls and instruments were started before the tunnel 
drive motor was started. The elevator was positioned so as to hold the 
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model on the mat, and tension was applied to the hold-down wire. When 
the tunnel airspeed reached approximately 100 mph, up elevator was 
applied until the model rose from the mat to approximately the center 
of the tunnel. The tunnel airspeed was then slowly increased and the 
elevator trim angle was varied manually to keep the model close to the 
t\mnel center horizontally. If the model deviated too far from the 
tunnel center laterally owing to drift in the roll-control gyroscope, 
the precession device was energized to return the model to the center. 
(The model rarely remained stationary in the tunnel so that trim control 
was applied only when it was evident that the model was oscillating about 
some location to the side of the tunnel center.) 

In general, the flights with the rigid wing model were terminated 
at the speeds at which the model motions were considered to he of such 
large amplitudes that the model was in danger of being destroyed. At 
such a time the model was landed after the tunnel airspeed was reduced 
to 100 raph. The gearing ratio of one of the autopilot components was 
then changed, and the model flight test repeated. This process continued 
on a trial basis until a set of autopilot gearing ratios was found which 
gave satisfactory flights at speeds up to 200 nqjh, the maximum speed of 
the tests. Some difficulty was encountered during this process in that 
the gearing ratios found to give good flights at low speeds were found 
to give unsatisfactory flights at high speeds. The converse trend was 
observed as well. However, compromise gearing ratios were found which 
permitted the model to be flown at both low and high speeds. 

For the model equipped with the flexible wings, the process was the 
same except that the maximum speed was the wing flutter speed. (The 
wing flutter speed (ref. l) was varied by moving weights in the external 
stores.) It should be mentioned that the flying characteristics of the 
rigid and flexible wing models were not identical (because of aero- 
elasticity) so that the gearing ratios used for the rigid wing model had 
to be readjusted for the flights with the flexible wing panels. 


CALCUIATIONS 


The eqxiations of motion for the towed model. Including the auto- 
pilot components, are presented in the appendix. As mentioned, the 
lateral equations of motion were obtained from reference 2. The aero- 
dynamic coefficients used in the investigation are listed in table III. 
The static derivatives were obtained using the experimental data of 
reference 5 and unpublished data, and the rotary derivatives were 
obtained using the procedures of references 6 to 8. Model mass param- 
eters (mass, moments of Inertia, center-of -gravity location, and incli- 
nation of the principal axis) were determined experimentally by swinging 
and weighing the model. 
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Lateral and longitudinal dynamic stability studies were made on an 
analog computer to obtain time histories of the model motions. An initial 
condition was used to start the motions. For the lateral motions, a 
sidewise displacement in y was used as an initial condition, whereas 
for the longitudinal motions, a vertical displacement z was used. All 
initial velocities, accelerations, and other displacements were ass\Jined 
to be zero. 

The roots of the lateral equations of motion were also obtained on 
an electronic calculator to define the periods and damping of the model- 
autopilot combination better. 


RESULTS AM) DISCUSSION 


Systematic calculations of the period and cycles to dang) to half 
ang)litude have been made to show the effects of the autopilot gearing 
ratios. The results of the analytic stijdy, for the lateral modes with 
three degrees of freedom of the rigid wing model at li)-5 ft/sec, are shown 
in figure 6. Three modes of motion are evident, and they are designated, 
for convenience, as the long-, middle-, and short-period motions. The 
effects of changes in K^i, Kp, and on the calculated periods and 

cycles to damp to half ang>litude axe shown. Some experimental periods 
are shown for cong)arison. For the case of K^, the curves are shown 
dotted in the range 1/4 < Yip < 5/4 since no points were calculated 
between those values and the dotted line represents the probable curve. 

In an attempt to approximate some of the model characteristics, a 
one-degree-of -freedom (in roll) analysis was made; the results of which 
are shown in figure 6. It can be seen that the middle period is fairly 
well predicted (for Yp between about 1/2 to 2), but the daflig)ing infor- 
mation shows a total lack of agreement with the three-degree-of -freedom 
analysis . 

lypical time histories of the rigid wing model motion are shown in 
figure 7* A power spectrum density analysis of the angle of roll 0 
was conducted in an attempt to find the natural frequencies. Since the 
motions were influenced by the tiinnel excitations, the natural frequencies 
were difficult to isolate from the forced frequencies. Figure 8 shows a 
typical power-spectrum-density analysis to illustrate the difficulties 
involved. However, the frequencies which were selected as being the 
natural frequencies are presented in table IV. The frequencies were 
chosen on the basts of this analysis and a judicioiis examination of the 
time histories to select aaiy relatively imdtsturbed portion of the data. 
The periods of the natural frequencies of both the lateral and longitu- 
dinal degrees of freedom are shown (periods for the longitudinal modes 
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were obtained from the time histories exclusively) . It can be seen that 
no general statement as to the reliability of the calculations can be 
formulated since the agreement varies from good to bad, and it is not 
clear whether the lack of correlation stems from inadequacies of the 
theoretical approach, inaccuracies in determining the stability deriva- 
tives, or from the inability to isolate the natural from the forced 
modes of motion. 

Figure 9 shows the motions of the flexible wing model \inder the 
influence of the tunnel disturbances. It can be seen that, in general, 
the an 5 )lltudes of the motions are somewhat greater than when the model 
is equipped with the rigid wing panels. 

As previously mentioned, the model apparently was continuously 
excited by imknown tunnel airstream disturbances during the flints, so 
that quantitative determination of the model damping characteristics 
was precluded. However, on the basis of visual observation, the oper- 
ators evaluated the effects of autopilot gearing ratio changes on the 
model damping. In general, these deductions confirmed the calculated 
effects of K^t and on the damping shown in figure 6. However, 

the calculated indications that changes in (above « 1') did 

not materially affect the cycles to damp were not observed. The operators 
felt that increases in improved the long-period damping. 

To illustrate the records obtained by an analog computer using the 
stability equations in the appendix, figure 10 shows the calculated time 
histories of the rigid wing model following an initial disturbance in 
y and z. The use of this initial condition was necessitated by the 
lack of knowledge of the actual forcing function, as previously discussed. 
However, the studies were of use for rapid determination of the trends of 
gearing ratios for stable flints. 


CONCLUDING REMARKS 


Exploratory wind-txmnel tests of a towed airplane model equipped 
with an autopilot demonstrated the ability of the autopilot system to 
control the model at speeds up to 200 n^jh. 

In seeking to make a comparison between the resiilts of the theoret- 
ical and experimental investigations, it must be recognized that the 
natioral frequencies may be obscxired by the forcing frequencies of the 
tunnel airstream disturbances, so that any differences in the results 
may be a result of the experimental data, as well as inaccuracies in the 
theory or the derivatives. Visual observation of the experimental trends 
qxxalitatively confirmed the calciilated variation of damping as the 
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autopilot gains K^t and were changed. However, the calculated 

effects of chajiges in were not visually confirmed. 

The calculations predicted three oscillations in the lateral degrees 
of freedom and two oscillations in the longitudinal degrees of freedom. 
These results were confirmed experimentally. 


langley Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
langley Field, Va., December l4, 195^* 


NACA RM L5^^L25 


15 


APPENDIX 

EQUATIONS OF MOTION OF A TOWED AIRPLANE 


The equations of motion of the towed airplane have been determined 
with respect to the stability axes (fig. l). In order to simplify the 
equations and to facilitate the mathematical solutions , the following 
assumptions are made: 

(1) All motions are perturbations from the equilibrium condition; 
thus, squares, higher-order terms, and products of the perturbations may 
be neglected. All angles are assumed small. 

(2) The airplane is initially flying with no roll or yaw and is in 
unaccelerated flight. 

(3) The basic aerodynamic parameters are constant during any motion. 

(4) The flow is quasi-steady; that is, any change in attitude pro- 
duces an instantaneous change in forces. 

(5) The towline is always straight. 

( 6 ) The towline is weightless and infinitely flexible and has no 
effect on the airplane other than that caused by direct forces and moments 
resulting from the towline tension. 

( 7 ) The lateral and longitudinal motions for the towed airplane are 
Independent and may be considered separately, as for the free airplane. 


Lateral Equations of Motion 

The lateral equations of motion for a towed airplane were derived 
in reference 2 and they are repeated here, in sli^tly different form, for 
completeness. 
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- CYp + 


_D^\ 

ZDb/ 


+ 



CYt 


- ^ I>b - Cw + Cd^ 0 = CYg 6r 




WzV - ^ D, . C, ^(l . ^ 


^ + 


-^b%z“b^ - ^ Eb + Cd ^)i* • * C„5^6a 


■ (Al) 




■^h^XTPb^ - ^ '*' f") ^ 


ZDb 


^ + 


^ Db + Cn^j0= Cls^Sr + Cjg^ 


The detenninant of the characteristic lateral-stability equations 
yields a ninth-order equation^ with three zero roots. After the zero 
roots are dropped^ the equation is of the form 

Db^ + ADb^ + BDb^ + CDb^ + EDb^ + FT»b + G = 0 (A2) 

For a detailed derivation of the equations, see reference 2. 

In addition to the airplane-stability equations, the following equa- 
tions were used to represent the various autopilot components: 


Roll control.- 


Sa = -y 


(A5) 
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Lateral au'toin&'blc position c on'fcrol plus yav dan^6 ^» 


&r = 


or. In operator form. 


6r = K^' 


1 + 


f) 


, (P + b 

^ I 


+ K^ir ^ 


(Alv) 


where the towline angle ilr' is approximated by 


,|r' 




P 4- i 

I>b 


for the frequency range of interest for these computations , the represen- 
Sl^cf the roil ^trol and yaw damper by ideal egjatlona 

is permissible since the response characteristics are linear in the low 

frequency range. 

The roll-control gyro and the yaw-damper rate -gyro were, of coi^se, 
fixed to the model. The quantities to which they responded were actually 

the roll about the body axis and the yawing L^ver, 

rather than roll and yaw about the corresponding stability axis, ^o , 

SS thTLSe of attack was Bmll, It was felt that It would 
the equations unnecessarily to transfer these Tihntltles 
axis system for this investigation. It was assumed 

autopilot components responded to the quantities about the stability axes. 


Longitudinal Equations of Motion 

The basic longitudinal equations of motion for a 
rtfsffrees of freedom were modified to Include the effect of the tow 
line forces and moments. No attempt has been made to analyze the ^ cs 

of other than to include the effect of the forces on the 

aiirlsine . 

Figure 11 shows a schematic representation of the towed ai^lane in 
a distiSJd corZTlln and will be used as an aid in developing the eq^- 
tiSs S Stion. It should be noted that the total instantaneous angle 

of attack 
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a' = Oq + a 

anH that the total instantaneous angle of pitch 

a = 0o + e 


and finally that 

= 0Q 

since the airplane is assumed to be in level flight in the undisturbed 
conditions. It should be noted that with the axes system being used 
(stability axes), the X-axis is fixed to one given position (the direc- 
tion of the undisturbed airstream) to the body during any disturbed 
motion. 


The longitudinal equations of motion for a towed airplane may be 
written as follows: 


Dc^e = CmD^QDcG + (CvoQ^aPc + Cmojo- + + 

pSc 

— ~ ~ _ Cx,p q'\d 00 - Cx^a, - RC^u - Cxjg &e + 

pSc \pSc =7 ® 

(cze)ie + (cZu)t“ * 

^ D^u = -Cjfi - (cDj^ - Ci)a + 2 Cbu + (%) ^» + 


(a5) 
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where the coefficients with the subscript T represent the forces and 
moments due to the towline tension. It is these quantities which m\ast 
be determined. 

To simplify the writing of the equations eind to transform them into 
familiar notation, let 


^^c = 


m 

pSc 




mc^ 


^mn 




_ ^ 

"®Dca 2 


2 


Substituting these values into equations (A5) leads to: 
2hcKy2d^20 ^ ^ + 

CmSgBe + (cmg) ^9 + (‘^)t^ (^z)t^ 


2p(.Dca = 1^2 " *"L8e^e + 

(%)t® * (%)t“ * (%)t“ 

2dcI^cU = -Cl0 - - Cl) a - 2Cpu + (Cxq) ^9 + 

(^Xa)T°^+ 
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The towline terms 
tute in equations (a6). 


Cx> ®ust be determined to substi- 

Referring to figure 11, it can be seen that 


A = sin-1 ^ 

L 


= sin-lA - 


or, with the sines replaced by the angles. 


A « 1 - - 


r “ ! - t(®c, e) 


(AT) 


Now the airplane was assinned to have initially been in level flight so 
that 


^ = A - 0 

^ • T t)® 

which is the angle between the towline and the X-axis. 

Consider the towline force along the X-axis: 

Tx = T cos I 

Within the approximation of small angles, this can be replaced by 

Tx « T (A9) 

Along the Z-axis^ the towline forces are 

T2, = -T sin I 

or again the approximation 


Tz - -T^ 


(AlO) 
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In normal level flight with the assumption of small disturbances, 
it can be assumed that 


T « D = Ci)Sq’ = Cd a'Sq' 

CL 


Equations (A9) and (AlO) now become 


Tk = CD^cc'Sq' 


(All) 


Tz = -CDa,°'’Sq'^ = -Cn^a'Sq’ 


z 

I 




(A12) 


Since the towline is attached to the nose of the airplane, at x', away 
from the center of gravity, the towline forces produce a moment about 
the Y-axis. This moment can be expressed as 


Mqi = -Tyx' sin 0Q - T^x' cos 6 q 


Mr -Txx'0o - T^x' 


(A13) 


Substituting equations (All) and (A12) into (A13) 


Mt ' (cDa“’ I + ='■ (cDa“’ | SV^) f ‘ T t) 


Mr = -x’ I SCDo,a'V^ 




(A14) 


The towline longitudinal force and moment derivatives can he found 
by differentiating equations (All), (A12), and (Aik) with respect to the 
desired variable. 
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Consider the X-forces: 
In coefficient form 



Similarly, 

Tx- = ^ = |-(ct) a’ £ V^s) = 0 

^0 Se Se V ®a> 2 / 

then 



0 


and the force due to velocity is 

’^V - ^ - 5?(ci>a“' I - Cb^»'pVS - CdpVS 

which in coefficient form is 



or 


K)o, = = 2Ci, 


(A15) 


(ai6) 


/ 


(AIT) 


(Al8) 
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then 




0 


(A19) 


It is recognized that for towed-model testing in a wind tunnel, V 
co\ild be neglected as a variable. For possible application to towed- 
gliders problems, V is retained as a variable at this time. 

Similarly, for the Z-forces: 








II 

-CDa“' 1 

T-T^o- 

CD 


= -CD, I V2S 





(A20) 


The use of this coefficient would, of course, resTilt in nonlinear 

equations since the coefficient is dependent upon both z and 0. Thus, 
to maintain the assumption of linear equations, let 




(A21) 


It shotild be noted that the z and 0 terms can be larger than the 
0Q term. 


2 h 
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Similarly, 


TZq 



Ct) a' 


then 


also, 



TZy 


^Tz ^ 
by by 



p 




-CD^a' P V2S 

CD 

o 

1 

- t ) 

0 

' “ -QDa'^'pVS 







then. 



2CD^a' 


V 




2 Cd 


Z 


V z 



or 


N)t ' ''K)t 



As before, this must be simplified to 

K)t = T 


I * t) 


(A22) 



(A25) 
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The force 


dT, 

= ^ = 


'•Z S 


1 

V2s 

1 

1 

o 

1 

0 



1— — J 



= -Cp a' 


I V2s 


= -CD I V2S i 


reduces to coefficient form as 



(A2l^) 


The pitching-moment derivatives due to the towline can be found from 
equation (A14) : 


MTa = 


^ 5a 


I v2sCD^a'x' 


0n + 91 + 


I 



z 


I 


or 



As before, neglecting the terms which will be products of small disturb- 
ances, so that linearity will be preserved. 



-Cd„|- eo 1 + 



(A25) 
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Similarly 


Mm - = i- 

do ^0 


- I V^SCDo^a'x' 


( 6 o e) (i + ^) - f 


I v2sCD„a'x' (l + . . I v%Cbx' (l + 


then 


Wt ' * t) 


(A26) 


also^ 


dMrr 

6 v 




- I 


9 o + 9 


-T 


= -pVSCc^a'x' 


f 9 o + ej[i + ^ -f 


r 


= -pVSCdx' 


-f 


then, 


N)t ■ ■ 


2 Cp 


'®o + 11 + T - f 


or 


= -2Cd f K®o 


J 
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Again neglecting potential products of small disturbances. 


(Cmu)^ = -2 Cd ^(l + y)6c 


(A27) 


when 


= 


^ . a_ 

Sz bz 


I v2sCD^a'x' 




= I v2Sx'CD^a' i 


the coefficient form is 


(<v) 


= Cn a’ 1 

z7t cl ° c 2 


(A28) 


Now inserting equations (A16) to (A19) and (A21) to (A2.Q) into eqtia- 
tions (a 6), the following expressions are obtained: 


2n Ky2Dc20 




Dc - ^ Cp 1 + ^ 


0 + 


D 4- C 


^ Cpa®o (•■■*■ t) ^ T ^ ~) 


u + 




2HcDc“ = 




2nc - + Cp 


e + Lci^ + Cp^ 2^ Q\a + 


^-2Cl + ^ 9qCp^u - ^ z - Cpggbe 

2hcDcU = -Cpe + (Ci, - Cp^ + Cp^^a + (2Cp - 2Cp)u 


(A29) 
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The last equation in equations (A29) can obviously be reduced so equa- 
tions (A29) become 








0 + 


Sx 


+ 


Sx 


f ^ T 


2L1 . n 5 


a + 


% - 2 2L.CD0O 1 + 


t) 


u + 




.2^Xo - 


Cl, 


2 ,»c ^ 1 . Y 


® ■*■ ("CLa + Cd^ — 0oja + 


2x' 


D 


-2Cl + 6oCd)^ “ ~ ^ ” Cjjg 


2|-^c®c^ = “Cl( 6 - a) 


' (A30) 


Letting 


= -V J(a - a 


')dt 


= -V 0 


^o “ ^ J ~ ^ 


)dt 


and remembering that the airplane was assumed to be in level flight, 
so that 0^0 = 00 ^ leads to 
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Then 


or in operator form 


-V / (do - ao)dt = 0 


= -V ~ 0')dt 


z = 


c(9 - g) 

D. 


Inserting this into equations (A50) gives 


2hcKy2Dc^9 = 




9 + 




-'a 


(1 + ^ + 2L1 2d 

a + 

\ ij I Dc 



s.- 




u + CniQgSe 


(A51) 


2 ^^Dca = 




(■ 


‘^®c + Cd 1 + ^ 


Cd 


c'\ , Cpc I 

I J I 


9 + 


2x' 


-Ci^ + Cd — 0o - Y p-la. + I-2 Cl + — eoC^lu - 


2ncDcU = -Cl (9 - a) 
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The characteristic longitudinal-stability eq^tion formed by the deter- 
minant of equations (A3l) is a sixth-order equation, with two zero roots. 
The fourth-order equation which remains is of the form 



+ HDc^ 


2 

ir>c + 


JD^ 


+ K = 


0 


As mentioned previously, the equations of motion can be simplified 
by considering the forward velocity to be a constant. With this assump- 
tion, equations (A 31 ) become 


2^XcKY2Dc20 



^ ^Da0o 





a + C 


“6. 


2ticDcO. 




+ 




CLse^e 


(A52) 


The autopilot component which was applicable to the longitudinal degrees 
of freedom was the elevator component of the automatic position control 
and the pitch damper. The equation used to represent these con 5 >onents is 


6e = Kq-AB* + K-0 


(A53) 


The towline angle in the vertical plane 0’ is 
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The Increment in elevator deflection due to the motions is proportional 
to the change in this angle due to the motion, or 

A0' = 0* - 0Q 


Therefore , 

A9' - -f (A3M 

and the autopilot -component equation becomes 

i)e - Kg, Kje 



In operator form, then 




(A35) 


It can be seen that the use of this autopilot-component equation does 
not alter the form of the fourth-order characteristic equation. 
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TABLE I 

MODEL CHARACTERISTICS 


Weight (mcxiel + tanks) , lb 


94.2 


Moments of inertia: 

I^, slug-ft^ 5 -774 

IYq i slug-f t^ 6 . 375 

IZq^ slug-ft^ 8.160 

no, ^ 

Center -of- gravity location measured from nose, ft (21 percent M.A.C.) 2.97 


Wing: 

Aspect ratio 

Sveepback of quarter-chord line, deg 

Taper ratio 

Area, sq ft 

Span, ft 

Mean aerodynamic chord, ft 

Airfoil section (normalHo c/4). . . 

Wing incidence, deg 

Dihedral, deg 

Wing tvd.st, deg 


. . . 5.45 

. . . 4o 
. . . 0.578 
. . . 9-02 
. . . 5.587 
. . . 1.673 

. NACA 64A010 

. . . 1.5 
. . . -3.5 
. . . 0 


Horizontal tail: 

Aspect ratio 3-59 

Sweepback of the quarter-chord line, deg 40 

Taper ratio 1*0 

Area, sq ft 1.55 

Span, ft 

Chord, ft 0.667 

Airfoil section normal to quarter-chord line NACA 64A009 

Tvfist, deg 0 

Dihedral, deg 0 

Tail length, ft 3.267 


Elevator : 

Hinge line, percent of stabilizer chord 70 

Chord, ft 0.200 

Inboard elevator span, ft O.518 

Outboard elevator span, ft O.518 

Elevator area (total), sq ft 0.207 

Vertical stabilizer: 

Aspect ratio 1.45 

Sweepback of quarter-chord line, deg 41.3 

Taper ratio 0.467 

Area, sq ft 1.2p 

Spam, ft 1.263 

Mean aerodynamic chord, ft O.927 

Airfoil section normal to quarter-chord line NACA 64A011 

Tail length, ft 2.928 


Rudder : 

Chord (parallel to fuselage center line), ft 0.228 

Span, upper segment, ft 0.55° 

Span, lover segment, ft O.55o 

Aileron: 

Hinge line, percent of wing chord perpendicular to quarter -chord line 75 

Aileron chord Inboard (parallel to center line) , ft 0.430 

Aileron chord outboard (parallel to center line), ft O.360 

Span (15.5 percent to 52 i>ercent b/2) , ft I.IO3 


Fuselage: 

Length, ft ... . 
Frontal area, sq ft 
Fineness^ ratio . . 


6.403 

0.5397 

7.75 
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TABLE II 

AUTOPILOT CHARACTERISTICS 


Roll Autopilot: 

Wei^t (total including relays), lb 10 

Gyro: 

Voltage, d-c 24 to 26 

Speed, rpm 15,500 

Torque motor: 

Voltage, d-c 24 to 26 

Stalling torque, In- lb 2.98 

Stalling current, amp 15* 7 

Gearing (motor to cam) 44.8:1 

Cam: 

Angle, deg 28 

Efficiency (approx.), percent 75 

Gearing (output arm to gimbal) 0.555:1 

Trim coil: 

Voltage, d-c 24 to 26 

Current at 24 V d-c, ma 290 

Torque on inner gimbal, in-oz 7 

Precession rate of outer gimbal (gyro speed, 11,700 rpm) 

(no load), deg/sec 1.5 


Damper Autopilot: 


Rotor : 

Speed, rpm - 

At 1 Ib/sq in 10,000 

At 7 1^/sq, in. . . • 90,000 

Natural frequency, cps - 

At 85,200 rpm 60 

At 20,000 rpm 155 

Percent critical damping - 

At 85,200 rpm 0.09 

At 20,000 rpm 0.21 


Pickoff (flapper-venturi nozzles) : 


Spacing 0.010” to 0.011" plus flapper 

Pressure, Ib/sq in 7 to 8 


Servo (ACS miniature pneumatic) : 

Stroke, in ±1/4 

Piston area, sq in. - 

Front 0.6012 

Rear 0.57^ 

Feedback spring (8.5 turns of 0.026 in. music wire), 

in. O.D 0.542 

Force output at 20 Ib/sq in. , lb 8 

Natural frequency, radians/sec 40 to 25 

Damping 0.6 to 0.8 


MCA RM L54L25 


55 


TABLE III 

AERODYMMIC COEFFIdENTS 


V, n?ih 

100 

200 


0.417 

0.104 


0.045 

0.025 

Cl<x 

5-555 

5.610 

■ 

0.264 

0.066 



-0.659 

-0. 688 



-0.741 

-0.701 



0.218 

0.197 



-0.095 

-0.092 



0.551 

0.582 

Cnr 

-0.555 

-0.552 

Clr 

0.119 

0.107 

CYp 

0.550 

0.111 

Cnp 

- 0 .o 64 

-0.026 

c, 

ip 

-0.209 

-0.294 



-4.492 

- 4.584 


-2.515 

-1.689 



0.060 

0.060 

^n&a 

- 0.024 

- 0.024 

^Ibr 

0.006 

0.006 

^“&r 

-0. 060 

-0.060 

CY&J. 

0.050 

0.050 



-0.555 

-0.549 



0.559 

0.168 
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TABLE IV 

COMPARISON OF EXPERIMENTAL AND CALCULATED PERIODS OF MOTION 
All periods are in seconds; velocity is in feet per second 


Wing 

V 

K^. 

K 0 


% 

^0 

Lateral Motion 

Longitudinal Motion 

Short 

Middle 

Long 

Short 

Long 

^exp. 

^calc . 

P 

■^exp. 

^calc. 

P 

exp. 

p 

calc . 

p 

exp. 

p 

^ calc . 

P 

■^exp. 

p 

calc . 

Rigid 

145 

1 

4.0 

0.48 

2.6 

0 

0.59 

0.73 

0.98 

1.15 

6.42 

8.67 

0.60 

0.58 

2.25 

2.16 

Rigid 

296 

.50 

5-5 

.50 

2.2 

.60 

.47 

.45 

1.55 

1.09 

4.50 

4.82 

.46 

.23 

1.62 

1.54 

Flexible 

162 

1 

5.2 

.54 

2.0 

.50 

.55 

— 

1.10 

— 

5.82 

— 

.92 

— 

2.10 

— 

Flexible 

280 

.50 

4.0 

.59 

2.0 

.50 

•59 

— 

— 

— 

2.74 

— 

.51 

— 

1.95 

— 

Rigid 

145 

.75 

5.6 

.50 

— 

— 

.72 

.59 

1.15 

1.10 

5.70 

9.45 

— 

— 

— 

— 

Rigid 

145 

.75 

4.0 

.50 

— 

— 

.69 

.68 

1.04 

1.11 , 

7.40 

9.23 

— 

— 

— 

— 
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A4 





Figiire 1.- The stability system of axes. Arrows indicate positive direc- 
tions of moments^ forces, and angles. This system of axes is defined 
as an orthogonal system having the origin at the center of gravity and 
in which the Z-axis is in the plane of symmetry and perpendicular to 
the relative wind in the undisturbed flight condition, the X-axis Is in 
the plane of symmetry and perpendicular to the Z-axls, and the Y-axis 
is perpendicular to the plane of symmetry. During disturbed flight, 
the axes remain fixed in the airplane. 



Figure 2.- Principal dimensions of the model. All dimensions in inches 

unless otherwise noted. 
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Figure 5*- Photograph of the roll-control autopilot component used in 

the towed model. 
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Automatic Position Control 



Pitch or Yaw Damper 


Tunnel vertical center line 




Down aileron 


Roll Control 

Figure 4.- Schematic representations of the con 5 )onents of the autopilot 

system used in the towed model. 




— tow- cable 



, flutter observer 



Figure Schematic representation of test arrangement in the Langley 

19 - foot pressure tunnel for the towed model tests. ^ 
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Figure 6.- Calculated lateral stability characteristics of the rigid 
wing model at V = l45 ft/sec, with varying autopilot gearing ratios. 
Symbols denote experimental periods. 
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(a) V = ft/sec> K^' = Ij = k.O} = 0A8> Kg . = 2.6; Kq* = 0. 
Figure 7 .- Experimental time history of the rigid wing model. 
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(b) ^ ^ 29 (> ft/sec; K^. = O.5O; = 5 . 5 > = 0.50; Kg. =2.2 


; Kq = 0.60. 


Figure 7 *- Concluded. 


Amplitude ratio of roll angle, db 



0 I 2 3 4 5 

Freguency, cycles /sec 


(a) V = ft/sec> = l.Oj = Ij-.O; 
= 0A8j Kq • = 2.6; K^- = 0. 



1 


0 1 2 3 4 5 

Frequency, cydes/sec 

(b) V = 296 ft/sec; K^. = 0-5; = 3-5j 

= O.3 O; Kqi = 2.2 ; Kq* = O. 60. 


Figure 8.- IJyplcal power- spectrum-density analyses of the angle of roll 0 of the rigid wing model. 
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(a) V = 162 ft /sec 


^,deg 


^,deg o 



S2,deg 



> ^ = 5-2} - 0.54; Kg. = 2.0} = 0.50. 


Figure 9.- Experimental time history of the flexible wing model. 
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(b) V = 280 ft/sec; K^. = O.5O; = 4.0; = 0.39; Kg, = 2.0; Kg = 0.50. 






Figure 9*- Concluded. 
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(a) V = 146 ft /sec; K^< 


= 1; = 4.0; = 0.48; Kq . = 2.6; Kg = 0. 


Figure 10.- Calc\ilated time history of the rigid wing model. 
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(b) V = 295 ft/sec; K^' = 0.50> = 5.5> = 0*30i Kq • = 2.2; Kq' = 0.60. 


Figure 10.- Concluded. 
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Figure 11.- Schematic representation of the model. Towline combination 
used in developing the longitudinal equations of motion. 
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